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We generalize Maxwell equations which describe the vacuum of quantum electrodynamics into the 
quantum form. This nontraditional approach is different from the widely used theory — -Quantum 
Electrodynamics. From another viewpoint, it could be a direction for interpreting quantum theories 
properly. 
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I. INTRODUCTION 



ll.lO.Ef 



Since quantum electrodynamics is quite mature, it is 
not widely known that Maxwell equations are first quan- 
tized quantum equations of a photon after the Planck 
constant h was multiplied to it 
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In fact, there has 
been an almost generally accepted way to change free 
noninteracting Maxwell equations into the form of Dirac 
equation The main purpose of this letter is to study 
Maxwell equations with self-interaction of a photon. 

Much attention was paid to study spinor formulation 
of Maxwell equations after Dirac [21 find the relativistic 
equation for a particle with spin-1/2. Landau and Peierls 
photon wave functions [3] could ensure that the modulus 
squared of wave function be the probability density of a 
photon, but this approach has some weaknesses as Pauli 
Pi had pointed out. And Laporte and Uhlenbeck Q gave 
Maxwell equations' spinor formulation 

ffii = 2(*2 + i*i), 

522 - 2(^2 

9i2 = 52i = 2*3, 

Qik is a spinor of second rank and ipj{ipj = Ej + iBj) 
is a complex 3-vector wave function whose real and vir- 
tual parts are electric and magnetic fields. The most 
crucial step was made by Silberstein @ who proved 
the possibility to describe Maxwell equations in terms 
of a 3-vector entities. However, it is Oppenheimer Q 
who noticed quantum nature of Maxwell equations for 
the first time(there is other opinion about this [1]). 
And Good 3 gave more clear interpretation about 
the similarity between Maxwell equations and Dirac 
equation, the new Maxwell equations transformation 
properties and conservation theorems. Following the 
idea of Riemann-Silberstein-Majorana-Oppenheimer ap- 
proach, Bialynicki-Birula pi]|4l2j gave the more detailed 
consideration about single photon quantum mechanics in 
resent years. 

In addition, Gersten p"3l - [l5| finds an alternative ap- 
proach to gain the Dirac form of the Maxwell equations, 
which starts from the first principle and the derivation 



process is just like Dirac equation. If the momentum p 
is written in form of 3- vectors column 



P = 



= [Pl P2 Pz] 



Decomposing the left side of the relationship between 
energy and momentum of electrodynamic field after mul- 
tiply a 3 X 3 imit matrix I, we obtain 
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E 



p^s 



E 



p^s 



PP 



0, 



where s = [si, S2, s^Y is a spin-1 vector matrix with three 
components 
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Then the photon equation is obtained, 



K - P'p) = f-I - P's) (-1 + p^s ) V- - ppV = 0,(2) 



where -0 is a 3 components wave function. At last, Eq. 
([2]) holds for the following two expressions are satisfied 
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c 



p^s 



0. 



(3) 
(4) 



If the quantum operator substitutions and the wave func- 
tion substitution are made as follows 



E ^ ifi^,V> =^ -ifiSJ, (5) 
V' = E-iB, (6) 

the classical Maxwell equations could be obtained 



-cs • p^/;. 



(7) 
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And the Gcrsten's approach return to Bialynicki-Birula's. 

In addition, there are many works ^16ri27j concerning 
the Dirac-like form of Maxwell equations. In this paper, 
we extend linear quantum theory of Maxwell equations to 
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nonlinear situation caused by the possibility of creating 
virtual particles in vacuum. Starting from Lagrange with 
an corrections in higher orders in E and B, we obtained 
the nonlinear Maxwell equations, and then try to convert 
it into the nonlinear Schrodinger equation with the same 
approach and wave function as linear situation. 



(15) 



Then put the I£ into Eq. ([14]) and Eq. dH]), nonlinear 
Maxwell equations is obtained 



II. NONLINEAR SCHRODINGER FORM OF 
NONLINEAR MAXWELL EQUATIONS 

The Lagrangian of Quantum electrodynamics in vac- 
uum was given by Heisenberg and Eulerjll]. It can de- 
scribe the phenomenon of optical birefringence and ex- 
periments of measuring the vacuum birefringence. 

In the limiting case of the stationary and homogeneous 
electromagnetic field, the exact expression of ^ would be 

where F^^, is the electromagnetic field tensor, and -F'"^ is 
the dual-field tensor contracted by F^j^^ with completely 
antisymetric unit tensor (the Levi-Civita tensor) e^"^^ . 
F^i, has the relation with as follow 

i^M- = ^A'^^ - duA^, (9) 
A^ is vector potential in the electromagnetic field, so 



F„, 



■ -El -E2 

El B3 

E2 -B3 

.E3 B2 —Bi 



—E3 
-B2 
Bi 




■yS 
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B2 


B3I 


Bi 





-E3 


E2 


B2 


E3 





-El 


B3 


-E2 


El 






(10) 



(11) 



For notation here, Greek indices and S nm from 

to 3. 

And ki, k2 determine the magnitude of the nonlinear 
correction 



ki 



4a^ 



fc2 = 



7a' 

45to4 



(12) 



45m4 

2 -. 

a is fine-sructure constant: a = -^-r- « -r^. 
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The last two terms of Eq. ([8]) are the corrections in 
B^ — E' and B ■ E which are Lorentz invariants, so ^ 
is invariable in any frame of reference. ££ satisfies the 
following Lagrange's equation 



Ja: 



5^ 



■5{d^,A,)' 



= 0. 



(13) 



If the Lagrange's equation holds for any situation, both 
terms must be zero. 



5^ 



0, 



(14) 



dt'E' = V X B', 
dtB = -V X E, 

V • E' 0, 

V • B = 0, 



(16) 
(17) 
(18) 
(19) 



where 



E' = (1 - kix)'E - fcayB, 
B' = (1 - fcia;)B -I- /csJ/E, 
a; = -2(E2-B2), 
y = -4E • B. 



(20) 
(21) 
(22) 
(23) 



Constructed with E and B, E' and B' take over their 
positions. 

In order to get the Dirac form of nonlinear Maxwell 
equations, we combine Eq. (14) and Eq. (15) 



(24) 



Substituting expressions (20) and (21) into Eq. we 
separate E and B after partial differential operators dt 
and V act on E' and B' 



"E'' 




B' 




= V X 


-iE 



Mdt 



'El- 




'iEi 


E2 




lE2 


E3 


= N 


1E3 


iBi 


Bi 


1B2 




B2 


jB3_ 




B3. 



(25) 



M and N are 6x6 matrices involving E, B and space 
partial differential operators. See Matrix elements of M 
and N in Appendix. 

In order to retain only time partial differential operator 
in the left side, both sides of Eq. (|25l) are multiplied by 



-El- 




-iEi 


E2 




iE2 


E3 




iE3 


iBi 


Bi 


lB2 




B2 


jB3_ 




B3. 



(26) 
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then ih times both sides of Eq. ((261 



ihdt 



-El- 

E2 
E3 
iBi 
iB2 

lB3 
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1 

1 









iE2 
iEs 
Bi 
B2 
B3. 











-1 
-1 





'El- 




E2 




E3 




iBi 




iB2 




iB^ 



(27) 



If we assume 
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(28) 



and 



■El- 
E2 
E3 
iBi 
1B2 
^Bs, 



(29) 



then, 



ihdti/j = Hip. 



(30) 



Regard Eq. ([30]) as the quantum equation of Maxwell 
equations including virtual processes (a photon the trans- 
form into a electron-positron pairs and interact with an- 
other photons). 

We decompose H into two parts: Hamiltonian Hq 
of classical electromagnetic fields and Hamiltonian Hint 
concluding virtual particles 



HiP = {Ho + H,nM. 



(31) 



Hq is almost identical to Eq. (16), which exists when the 
wave function has 6 components, 
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Hq can be expressed with s and p 
Ho = 



(32) 



s-p 
s • p 



(33) 



Pk is momentum operater and Eq. (1) has given the ex- 
actly form of Sfc. Hint could also be written in a compact 
manner 



Hi, 



(34) 



where are 6x6 block matrices and consist of 3 x 3 
matrixes £k and Bk except zero elements 



—£k Bk 




(35) 



Elements of £k and Bk are given in appendix and they 
satisfy the following relations 



£k 



S^^(E,B)=-Sf(B,E) 



(36) 
(37) 



where fc run from 1 to 3. Eq. (37) means that, if 
we exchange E and B in the two elements which are 
symmetrical about the diagonal of T-Lint and reverse its 
sign, the two elements are equal. 

In order to gain the same wave function as Eq. (6), 
transformation is made 



ihdtU 



E 
iB 



UHqU-^U 



E 

iB 



UHntU-^U 



E 

iB 



where U is the transformation matrix 



V2 



1 1 

1 -1 



and U ^ is the inverse matrix of U. Then Eq. 
be written 



ihdt 



where "Ho and Hint is the new Hamiltonian 



"E + iB 


— Ho 


E + iB 


+ Hint 


■e + 


iB 


E iB 


E iB 


E 


iB 



Ho = UHnU-^ = 



s-p 
s-p 



Hint — -7: 



-£k + Bk -Ek-Bk 
^Ek + Bk -Ek-Bk 



■Pk- 



(38) 

(39) 
can 

(40) 

(41) 
(42) 



Obviously, the linear part of the Hamiltonian Ho, which 
describes the classical electrodynamics, is a Hermitian 
operator. However, Hint is not a Hermitian operator al- 
though the relation between its elements shown in Eq. 
(36) and Eq. (37) is enlightening. Based on the Hamil- 
tonian of Heisenberg and Euler which describes two pho- 
tons scatter off one another in empty space (virtual par- 
ticles are considered), we have got the Hamiltonian in- 
cluding Hint- And the Hamiltonian is apparently not 
Hermitian. Accordingly, the number of photons is not 
conservative in the colliding process of two photons'. 
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III. CONCLUSION 

In conclusion, starting from the Lagrangian ^ of 
quantum electrodynamics in vacuum, we give the cor- 
responding Maxwell equations. Then a new approach 
which is different from the tradional method to derive 
from the form of nonlinear Schrddinger equation. At 
last, we transform the wave function in order to obtain 
Riemann-Silberstein-Majorana-Oppenheimer wave vec- 



tor representation. 



Acknowledgments 

This work is supported in part by NSF of China 
(Grants No. 10975075), Program for New Century Ex- 
cellent Talents in University, and the Project-sponsored 
by SRF for ROCS, SEM. 



[1] W. Greiner, Relativistic Quantum Mechanics, (Springer- 

Verlag, Berlin, 1990), pl05. 
[2] P. Dirac, Proc. Roy. Soc. 117, 610 (1928). 
[3] L. D. Landau and R. Peierls, Zeitschrift fiir Physik A 

Hadrons and Nuclei, 62, 188 (1930). 
[4] W. Pauli, Prinzipien der Quantentheorie, Handbuch 

der Physik, Vol. 24, Springer, Berlin, 1933, (English 

translation: General Principles of Quantum Mechanics, 

Springer, Berlin, 1980). 
[5] O. Laporte, G. E. Uhlenbeck, Phys. Rev. 37, 1380 (1931). 
[6] L. Silberstein, Annalen der Physik, 327, 579 (1906). 
[7] J. R. Oppenheimer, Phys. Rev. 38, 725 (1931). 
[8] E. Mignani, E. Recami, M. Baldo, Lettere Al Nuovo Ci- 

mento, 11, 568 (1974). 
[9] R. H. Good, Phys. Rev. 105, 1914 (1957). 
[10] L Bialynicki-Birula, Acta physica Polonica. 86, 97 

(1994).' 

[11] L Bialynicki-Birula, Acta Physica Polonica B, 37, 935 
(2006). 

[12] L Bialynicki-Birula, Phys Rev A, 79, 032112 (2009). 
[13] A. Gersten, Ann. Fond. L. de Broglie, 21, 67 (1996). 
[14] A. Gersten, Foundations of Physics Letters, 12, 291 
(1999). 

[15] A. Gersten, Foundations of Physics, 31,1211 (2001) 
[16] T.Ohmura, Progress of Theoretical Physics, 16, 684 
(1956). 

[17] R. Cook, Phys. Rev. A. 25, 2164 (1982). 
[18] R. Cook, Phys. Rev. A. 26, 2754 (1982). 
[19] E. Giannetto, Lett. Nuovo Cim. 44,140 (1985). 
[20] H. Nuez, A. Salas Brito, A. Salas Brito, C. Vargas, Re- 
vista Mexicana de Fisica. 34, 636 (1988). 
[21] T. Inagaki, Phys. Rev. A. 49, 2839 (1994). 
[22] J. Sipe, Phys. Rev. A. 52, 1875 (1995). 
[23] S. Esposito, Found. Phys. 28, 231 (1998). 
[24] T. Ivezic, Found. Phys. 35, 1585 (2005). 
[25] T. Ivezic, Found. Phys. Lett. 18, 401 (2005). 
[26] T. Ivezic, EJTP, 10, 131 (2006). 

[27] A. Bogush, V. Red'kov, N. Tokarevskaya, G. Spix, 

larXiv:0905.0"26TV l (2009). 
[28] W. Greiner, Quantum Electrodynamics, (Springer- 

Verlag, Berlin, 2002), p426. 



Appendix A: Elements of £k and Bk 



Si ■■ 


Si' 


= 0, 




: 


= -8kiE2B3+8k2E3B2, 




: 


= 8kiE3B2 - 8k2E2B3, 




sr- 


= £f = -AkiEiB3 + Ak2E3Bi, 




sl' 


= £f =AkiEiB2-4k2E2Bi, 






= £f = 4((fci - k2){E2B2 - E3B3) 


S2 ■■ 


cU 
^2 


^ 8kiEiB3 - 8k2E3Bi, 




c22 
^2 


^ 0, 




c33 
^2 


= -8kiE3Bi - 8k2EiB3, 




cl2 
^2 


= £f = AkiE2B3 - Ak2E3B2, 




cl3 
^2 


= £f = 4(fci - k2){E3B3 - EiBi), 




1723 
^2 


= £f = -4kiE2Bi + 4k2EiB2. 


S3 : 


^3 


= -8kiEiB2+8k2E2Bi, 




c22 
^3 


= 8kiE2Bi ~ 8k2EiB2, 




C.33 
^3 


-0, 




^12 
^3 


- £f = 4(fci - fc2)(-BiBi - E2B2), 




<7l3 

^3 


= £f ^ ~4kiE3B2 + 4.k2E2B3, 




c23 
^2 


= £f = 4kiE3Bi - 4k2EiB3. 


Bi : 


Bl' 


= 0, 




Bf 


= ^ki-k2){B2B3- E2E3), 




Bf 


= A{ki - k2){E2E3 - B2B3) 




BP 


= -4kiEiE3 - Ak2BiB3, 




Bf 


= AkiBiB3 + Ak2E3Bi, 




B\^ 


^AkiEiE2+4k2BiB2, 




Bf 


= ~4kiBiB2 ~ 4k2EiE2, 




Bf 


= ikfEl + Bl) + Ak2{El + Bl), 




Bf 


= -AkfEl + Bf) - Ak2{El + Bl). 



(Al) 



(A2) 



(A3) 



(A4) 



0, 

-4{h - k2){EiE3 - B1B3) 
-AkiB2Bz - Ak2E2Ei, 
AkiE2E3 + Ak2B2B3, 
-4ki{El + Bl) - Ak2{El + Bl), 
Aki{El +Bl)+ 4k2{Ef + Bl), 
-AkiEiE2 - Ak2BiB2, 
4kiBiB2 + Ak2EiE2. 

-4{ki-k2){E^E2-B,B2), 
4(fci - k2){EiE2 - B1B2), 
0, 

AkiiE^ +Bi)+ 4k2{Ei +Bj), 
-Aki{El + Bf) - 4k2{Ef + Bl), 
AkiB2B3 + Ak2E2Es, 

-4kiE2E3 - 4k2B2B3, 

-AkiBiBs - 4k2EiE3, 
AkiEiE3 + Ak2BiB3. 



